Static Quark Potential and the Renormalized 
Anisotropy on Tadpole Improved Anisotropic 

Lattices 1 



Wei Liu a , Ying Chen b , Ming Gong a , Xin Li a , Chuan Liu' 

and Guozhan Meng a 

^School of Physics, Peking University 
Beijing, 100871, P. R. China 

h 'Institute of High Energy Physics 
Academia Sinica, P. 0. Box 918 
Beijing, 100039, P. R. China 



Abstract 

Static quark potential is studied using a tadpole improved gauge lattice action. The 
scale is set using the potential for a wide range of bare parameters. The renormalized 
anisotropy of the lattice is also measured. 
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1 Introduction 



It has become clear that anisotropic lattices and improved lattice actions 
greatly facilitate lattice QCD calculations involving heavy physical objects 
like the glueballs, one meson states with non-zero three momenta and multi- 
meson states with or without three momenta. It is also a good workplace for 
the study of hadrons with heavy quarks. In many of the lattice calculations, it 
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is useful to use finite size techniques in which lattices with fixed physical vol- 
umes are simulated. For the purpose of these calculations, one needs to know 
the correspondence between the lattice spacing in physical units and the bare 
gauge coupling (3. It is also important to check the renormalization effects of 
the anisotropy parameter for these actions. In this work we present our nu- 
merical studies on these issues using the tadpole improved gluonic action on 
asymmetric lattices: 

BTrPij 1 TrRij 1 TrRji 
9~&4 ~ 36 Z u e s ~ 36 £ u* _ 

± tjoTrP 0i I &JVjjo l m 
9 u 2 s 36 \ U 

where is the usual plaquette variable and Rij is the 2x1 Wilson loop on 
the lattice. The parameter u s , which we take to be the forth root of the av- 
erage spatial plaquette value, incorporates the usual tadpole improvement [1] 
and £ designates the bare aspect ratio of the anisotropic lattice. With the 
tadpole improvement in place, the bare anisotropy parameter £ suffers only 
small renormalization. Therefore, the physical, or renormalized anisotropy pa- 
rameter £, which can be determined by comparing physical quantities which 
depend on distances in both the spatial and the temporal directions, will be 
quite close to its bare value £o- The effect of this renormalization will be stud- 
ied in this paper for various values of f3 using the so-called side-way potential 
method. Using the pure gauge action (1), glueball and light hadron spectrum 
has been studied within the quenched approximation [2,3,4,5,6,7,8]. 

In this paper, we will study the static quark anti-quark potential using the 
pure gauge action given in Eq. (1). The static quark anti-quark potential is 
obtained by the measurement of the Wilson loops on the lattice. Then, we 
set the scale in physical unit using the Sommer scale r which is defined via 
the static quark potential. This establishes the relation between the lattice 
spacing and the bare gauge coupling. For all the (3 values we studied, we also 
investigate the renormalized anisotropy parameter £ using the side-way po- 
tential method. Similar studies have been carried out before [9,10]. Our study 
covers more parameter space which is used in anisotropic lattice simulations. 
The renormalized anisotropy is also calculated within perturbation theory to 
one-loop [11]. 

This paper is organized in the following manner. In Section 2, parameters of 
our simulations are given and the details of the Wilson loop measurements 
are presented. From these data, we extract the static quark potential and the 
lattice spacing is determined in terms of the physical scale r$. The results for 
all (3 are then interpolated using quadratic polynomials. These interpolation 
then offers a rather precise correspondence between (3 and r /a s over the whole 
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range of (3 that have been simulated. In section 3, we discuss the side- way 
potential measurements which yield the renormalized anisotropy for all our 
simulation points. In Section 4, we will conclude with some general remarks. 



2 Simulation Results for Static Potential 

In this section, we present our numerical results for the study of the static 
quark anti-quark potential for various gauge coupling j3. Our values of (3 range 
from 1.9 to 2.8 which roughly corresponds the lattice spacing a s in the range 
of 0.35fm to 0.12fm. Details of the simulation parameters are summarized in 
Table 1. 

Since we are interested in the spatial lattice spacing in physical unit, which is 
sensitive to the ultra-violet physics, we therefore only use small lattices of size 
8 3 x 40. The lattices are generated using the conventional Cabbibo-Mariani 
pseudo-heatbath algorithm with over-relaxation. In a compound sweep, we 
perform one pseudo-heatbath sweep with three over-relaxation sweeps. For 
each p, over a thousand of gauge field configurations are generated, with each 
configuration separated by 3 compound sweeps mentioned above. 

Using these gauge field configurations, Wilson loops are constructed and mea- 
sured which yield static quark anti-quark potential. It is known that, Wilson 
loops, especially large loops are extremely noisy objects. We thus utilize the 
conventional smearing techniques for the spatial gauge links. In such a process, 
each spatial gauge link is replaced by its original value plus a linear combina- 
tion of its nearest staples with a coefficient X s . Finally, each spatial gauge link 
is projected back into group 577(3). Two different sets of smearing parameters 
are chosen for small spatial distance R and large R. Typically the smearing 
parameter in our simulations lies in the range: X s = 0.1 — 0.4. Smearing can 
be done iteratively. We usually perform 4 or 6 sweeps of smearing before our 
construction of the Wilson loops. Another useful technique is to use thermally 
averaged temporal links. This also greatly reduces the statistical errors. 

The Wilson loops W(R,t) at a fixed spatial distance R and large temporal 
separation t is related to the static quark anti-quark potential V(R) by: 

W(R, t) ~ e - tv W = e~ iV{R) , t -> oo . (2) 

Note that here we have used dimensionless temporal separation i — t/at which 
assumes integral values. What is really measured is the dimensionless poten- 
tial: V = a t V(R). Within the range that we can reach in our simulation, we 
find that the potential can be well-represented by a linear term plus a Coulomb 
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term. Therefore, at various values of R, we fit the potential V(R) using the 
following form: 

V(R) = V + ^ + aR, (3) 



Or, in terms of the dimensionless potential: 

* ot 

V(R) = V + ^ + aR, (4) 



with Vo = <kVo, a = aa t /a s , and a = aa t a s . 

At this stage, we tried two ways of extracting the parameters a and a. In the 
first method, one first checks for plateau behaviors of ratios of Wilson loops 
between two neighboring t in the large temporal region. According to Eq. (2), 
the height of these plateaus will give us the estimate for the dimensionless 
potential V(R) for all values of R. Then, in the second step, one fits the 
resulting V(R) versus R using Eq. (4). The result of this fit then gives the 
optimal values for a and a. In the second method, one directly performs a 
correlated fit for the on-axis Wilson loops using the form: 



W(R,t) = Z(R)exp 



a 



-t IVo + ^ + aR 



(5) 



Note that here all the Z(R) parameters plus the three parameters: Vo, ot and 
a enters the fitting process. 

In Fig. 1 and Fig. 2, we show the fit of the static potential for two differ- 
ent values of (5. The dashed curves are the fits for the potential using the 
first method. For comparison purposes, potentials using the fitted parameters 
obtained with the second method are also plotted as the solid curves. The 
difference for the two methods is indistinguishable for large values of (3. In 
any case, we find that the two methods always yield compatible results for the 
fitted parameters. We list in Table 1 the fitted parameters a and a for all f3 
values together with the y 2 per degree of freedom of the fit. 

To set the scales of a lattice in physical unit, we use the so-called Sommer 
scale r defined as: 



2 dV(R) 



R dR 



= 1.65 . (6) 
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Fig. 1. The dimensionless potential V{R) is fitted versus R using Eq. (4) for (3 = 1.9. 
The dashed curve is the fit using the first method while the solid curve represents 
a direct fitting to the Wilson loops. 
Table 1 



Simulation parameters used in this work. All lattices are of size 8 3 x 40 and £o = 5. 



(3 


1.9 


2.0 


2.1 


2.2 


2.3 




0.328 


0.345 


0.361 


0.377 


0.394 


a 


-0.017(2) 


-0.021(2) 


-0.027(2) 


-0.030(1) 


0.040(1) 


<7 


0.157(1) 


0.1335(9) 


0.1123(7) 


0.0937(6) 


0.0733(6) 


X z /d-o.f 


1.45 


1.04 


1.58 


0.89 


0.48 


ro/a s 


1.415(12) 


1.522(12) 


1.644(16) 


1.787(15) 


1.991(18) 


A 


14(2) 


14(2) 


14(2) 


17(2) 


20(2) 


B 


-44(7) 


-44(7) 


-44(7) 


-58(9) 


-71(11) 


C 


39(8) 


39(8) 


39(8) 


54(9) 


69(12) 


P 


2.4 


2.5 


2.6 


2.7 


2.8 




0.409 


0.424 


0.437 


0.451 


0.463 


a 


-0.0389(9) 


-0.0407(6) 


-0.0410(3) 


-0.0599(5) 


-0.0582(5) 


<7 


0.0606(4) 


0.0482(3) 


0.0384(2) 


0.0282(2) 


0.0225(2) 


X z /d.o.f 


0.90 


1.94 


2.81 


1.88 


0.83 


ro/a s 


2.191(17) 


2.449(17) 


2.743(11) 


3.093(20) 


3.476(23) 


A 


23(2) 


31(3) 


34(3) 


34(3) 


34(3) 


B 


-88(11) 


-125(15) 


-142(16) 


-142(16) 


-142(16) 


C 


88(13) 


135(19) 


155(20) 


155(20) 


155(20) 



This implies that: 



r fa + 1.65 _ la + 1.65/f 
a s Y aa 2 s V o 
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Fig. 2. Same as Fig. 1, but for (5 = 2.6. 

where £ = a s /at is the (renormalized) anisotropy of the lattice. As we will 
verify, the renormalization effects for £ are quite small, therefore we can re- 
place £ in the above expression by £o i n most cases. Also in Table 1, we list 
the values of (ro/a s ) for various f3 values. 

In order to obtain a one-to-one correspondence between (3 and r /a s within 
the whole range we studied, we have also attempted to make interpolations 
within the whole range of j3 in our simulation. We use a five point quadratic 
fit in 1//3 around each point to interpolate points between the neighboring 
points of f3 that have been simulated: 

a s (3 (3 2 



The situation is shown in Fig. 3. The data points are the resulting values of 
ro/a s obtained from the static quark potential at various values of (3. The 
colored lines are the quadratic fitting curves around each simulated point. 
We see that these quadratic interpolations are smoothly connected with one 
another. Within each interval, different fitting curves agree within errors. For 
reference, we tabulated the fitted parameters A, B, C for each interpolation 
in Table 1. With these parameters, one can find the correspondence between 
(3 and the physical lattice spacing within the whole range of interest. 
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1 I I I I I I I I I I I I I I I I I I I I 

0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50 0.52 0.54 

1/p 

Fig. 3. Interpolations in the whole range of (3 that we simulated. Around each simu- 
lated point, we perform a five point quadratic interpolation around that point. These 
interpolations are shown as colored lines. It is seen that these lines interconnect each 
other well, showing that consistency of the interpolation. 

3 Determination of renormalized anisotropy 

We now move to the measurements of the anisotropy parameter £ = a s /a t . 
This is a physical quantity which can be measured by various methods. In this 
work, we use the side-way potential method. The strategy for this method is 
briefly described below. 

On anisotropic lattices, one can measure both the spatial-temporal Wilson 
loops and the spatial-spatial Wilson loops. Both of these objects contain in- 
formation about the static quark anti-quark potential as a function of the 
quark distances. Since there is one unique physical potential, by comparing 
the results from these two measurements, one infers the information about 
physical anisotropy. 

To be more precise, we assume that the spatial lattice axis are called x, y and 
z, the temporal axis is called t. We can measure a spatial-temporal Wilson 
loop with the spatial separation R = Rz lies in the z direction, with z being 
a unit vector in the z direction. Since we are on a lattice, it is obvious that 
R = Ra s , with R being positive integers. If we measure the spatial-temporal 
Wilson loop W(Ra s ,ta t ) and we investigate the large R behavior of the loop, 
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we have: 



W(Ra s ,ta t ) oc e - kasVt{tat) , R > 1 . (9) 

Therefore, at large values of R, we obtain the combination a s Vt(iat) for all 
values of t. Here we use the notation V t to indicate that it is a potential 
obtained from a spatial-temporal Wilson loop measurement. 

On the other hand, we can also measure the spatial-spatial Wilson loop 
W(Ra s , R'a s ) with R'a s indicating the distance between the quark and anti- 
quark separation in the xy plane. We thus have: 

W(Ra s , R'a t ) oc e -A*.Vi(A'a.) ( 10 ) 

Again, for large enough R, we obtain the potential a s V s (R'a s ). Here we use 
V s to indicate that it is the potential from a spatial-spatial Wilson loop mea- 
surement. Note, however, Vt and V s are in fact the same physical quantity. 
That is to say, if they evaluated at the same physical distance, they ought to 
be identical. Therefore, by measuring the spatial-spatial and spatial-temporal 
Wilson loops, we obtain two versions of the same potential. The renormalized 
anisotropy £ has to be such that: 

V(R'a s ) = vda,) , (11) 



holds for all values of R! . 

To determine the anisotropy £, we plot a s V(^a s ) at all values of t and inter- 
polate in the whole range of i. With this interpolation, a curve is obtained 
which gives the values of a s V(^a s ) versus t continuously within the range. In 

the same plot, we draw horizontal lines with height of a s V(R'a s ) at various 
values of R' . The crossing of these horizontal lines with the interpolation curve 
then offer the value of t such that Eq. (11) is satisfied. This in turn gives an 
estimate for the value of £. This situation is demonstrated in Fig. 4. The data 
points are the results of a s V(^a s ) at different values of t. The interpolation 
curve is also shown. Since our temporal lattice spacing is finer, usually linear 
interpolation is adequate. The horizontal dashed lines correspond to the val- 
ues of a s V(R'a s ) at various values of Rl . The crossing points then finally give 
the results for the renormalized anisotropy. 

The anisotropy obtained in this way depend weakly on the values of R! chosen 
to make the connection. We choose to quote values of £ obtained at RI — 2 
as our final estimates since they are most stable and have less errors. Finally, 
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(3=2.6 




time,t 



Fig. 4. Determination of physical anisotropy using the side-way potential method. 
This is the case for (3 = 2.6. 

we can determine the physical anisotropy £ for all values of (5 that have been 
simulated. In Fig. 5, we plot the ratio of the renormalized anisotropy to the 
bare anisotropy: r\ = £/£ versus (3. It is seen that the relative renormalization 
of the anisotropy parameter is very mild, less than a few percent, within the 
whole range of f3. This is mainly due to tadpole improvement. It is known that 
at intermediate values of bare gauge coupling, the renormalization can be as 
large as 30% without tadpole improvement. 

It is also interesting to compare our measured values with the result from 
perturbation theory. According to Ref. [11,12], the renormalized anisotropy 
can be well fitted with the formula: 



7/ 



/ nnnpp 0.0702 0.0399\ 2 

1 + 0.0955 -— -g- )gl 

\ Xo Xo / 



(12) 



where the "boosted" bare coupling g% and the boosted bare anisotropy xo are 
given by: 

13 = — — , Xo = £o— • 13 

% Ut 



In our simulation, we find u t — 1 and when plugged into the above formulae, 
qualitative agreement between our measured values and the perturbative re- 
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Fig. 5. The effects of renormalization for the anisotropy is shown. We plot i] = £/£o 
at various values of (3, ranging between 1.9 and 2.8. 



suits are found. However, our measured values of 77 are closer to unity than 
the predictions from perturbation theory. 



4 Conclusions 



In this paper, we present a numerical study of the static quark anti-quark 
potential using anisotropic lattice gauge action. Using the static quark poten- 
tial obtained from Wilson loops, we determine the physical scale of the lattice 
spacing for a range of [3 values. This will set up a direct correspondence be- 
tween the value of j3 and the physical lattice spacing. We also measure the 
renormalized anisotropy of the lattice using the side-way potential method. It 
is found that the renormalization effect for the anisotropy is small over the 
whole range of (3 that have been studied, typically below 5%. This is due to 
the use of tadpole improvement. The results obtained in this study will be 
useful for further applications of anisotropic lattices on other issues like the 
hadron-hadron scattering. 



-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 
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1.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50 0.52 
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